Dynamics of magnetic domain wall motion after nucleation: 
Dependence on the wall energy 
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The dynamics of magnetic domain wall motion in the FeNi layer of a FeNi/AbOa/Co trilayer has 
been investigated by a combination of x-ray magnetic circular dichroism, photoelectron emission 
microscopy, and a stroboscopic pump-probe technique. The nucleation of domains and subsequent 
expansion by domain wall motion in the FeNi layer during nanosecond-long magnetic field pulses 
was observed in the viscous regime up to the Walker limit field. We attribute an observed delay of 
domain expansion to the influence of the domain wall energy that acts against the domain expansion 
and that plays an important role when domains are small. 

PACS numbers: 75.60.Jk, 75.60.Ch, 75.70.Kw, 75.50.Bb 
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Many of the fascinating phenomena in thin fihn 
magnetism that appear promising for new apphcations 
in magnetic data storage involve the use of muhilay- 
ered magnetic systems, for example spin-valves 1] and 
magnetic-tunnel-junctions "S]. The study of the time de- 
pendence of the magnetic properties in multilayered sys- 
tems is of high actual interest with respect to the achiev- 
able magnetic response time of such devices, but also 
from a fundamental point of view. It has been observed 
that the magnetization reversal for hard/soft magnetic 
trilayers is dominated by domain wall motion in quasi- 
static conditions and by the nucleation of many small do- 
mains in fast dynamic conditions 0, Q • Such a difference 
in the reversal mechanism must obviously be accompa- 
nied by dramatic changes in the energetics, and magnetic 
domain walls will become much more influential at faster 
reversal rates. 

The role of domain wall energy on the magnetization 
reversal is a long-debated issue. It has been discussed 
in the seventies in connection with bubble domains in 
perpendicularly magnetized materials. Malozemoff and 
Slonczewski pointed out that the domain wall energy 
leads to a wall curvature pressure that tends to contract 
an isolated bubble domain, eventually collapsing the do- 
main when the opposing demagnetizing energy is over- 
come 0. Skomski et al., in a self-consistent magnetic- 
viscosity approach, pointed out a negative velocity for 
bubble domain expansion in small fields as consequence 
of the domain wall energy, thus confirming the domain 
collapse 0. More recently, Lyberatos and Ferre studied 
the rate of domain growth in an applied magnetic field 
by Monte Carlo simulations. The results showed a slow- 
ing down of domain wall motion and the appearance of a 
fractal shape of the domains due to pinning effects in the 
regime of thermally activated domain wall motion. How- 
ever, no influence of domain wall energy was found in the 



viscous motion regime, in which pinning processes do not 
hinder the domain wall propagation Experimentally, 
it has been shown that the influence of the domain wall 
energy can lead to a straightening of domain walls in the 
low velocity creep-like motion regime in the presence of 
an introduced defect 8], and to a selectively accelerated 
domain expansion of such wall motion processes leading 
to a merging of domains JJ. 

In this Letter we present an experimental investiga- 
tion of the influence of the domain wall energy on the 
magnetization reversal of the soft magnetic layer in mag- 
netic tunnel junction-like trilayers. We study the reversal 
of the soft layer away from the magnetization direction 
of the hard layer. In this configuration the reversal is 
nucleation-dominated reversal for high fields 0, , and 
a strong infiuence of the domain wall energy is exp ected. 
We use a stroboscopic pump-probe technique |lCI| to ob- 
serve in real space the repeated nucleation and expansion 
of domains and find that the infiuence of the domain wall 
energy leads to an apparent delay in domain nucleation, 
which can amount up to several nanoseconds. This can 
be well described by an effective specific domain wall 
energy exerting a negative domain wall pressure when 
domains are small, just after nucleation, but becoming 
negligible once domains exceed a certain size. 

In in-plane domains there is no stabilizing demagne- 
tizing energy as in perpendicular bubble domains, and 
the domain wall energy P^d of a domain with perime- 
ter length P in a film of thickness d can be described 
simply by an effective field Hyjaii = 1^/ {^IJ-oMsP) acting 
against the external field H^xt Here 7 is an effective 
domain wall energy averaged over the perimeter of the 
domain, and A is a geometric factor describing the ra- 
tio between area and squared perimeter of the domain 
(A = Ait for a circular domain). Only if H^aii is small 
enough so that the corresponding wall velocity can be 



neglected, domains are stable in static conditions. In our 
case both, the domain wall field and the effective field 
from the coupling between the two ferromagnetic layers 
in a trilayered structure, Hcopi, oppose the externally ap- 
plied field, so that the total effective field is 



Htnt — Hp, 



Hcopl 



7A 



(1) 



and thus depends on P. We consider as measurable 
quantity the speed of perimeter expansion, dP/dt, which 
in turn will be a function of Htot- We show that the 
resulting inhomogeneous non-linear differential equation 
for P can well describe the observed apparent delay in 
reversible domain nucleation using a fixed value for 7A, 
thus revealing the manifestation of the domain wall en- 
ergy in the reversal dynamics. 

The sample used in this study was a 4 nm Fe2oNi8o /2.6 
nm AI2O3/7 nm Co trilayer grown by RF sputtering on 
the step-bunched surface of a 6° miscut Si(lll) crystal. 
Underneath the Co layer, 3 nm of CoO has been de- 
posited to increase the coercivity of the Co layer. The 
FeNi layer has been capped by 3 nm of Al to prevent 
oxidation. Details of the thermal treatment leading to 
the creation of step-bunches are given in Ref. [13. Ac- 
cording to quasi-static longitudinal Kerr effect experi- 
ments, the sample features uniaxial anisotropy along the 
step-bunches. The coercivities are 1 mT for the FeNi 
and 4.3 mT for the Co layers, and a Neel type cou- 
pling field fioHcopi of 1.2 mT favors parallel alignment. 
Details of similar types of samples are shown in Refs. 

The experiment was performed using a combination of 
x-ray magnetic circular dichroism (XMCD), photoelec- 
tron emission microscopy (PEEM), and a stroboscopic 
pump-probe technique jlO| . Magnetic field pulses are pe- 
riodically applied at 312.5 kHz by a coil made of Cu foil 
wrapped around the sample. The magnetic response is 
stroboscopically probed by XMCD-PEEM with the same 
frequency at a fixed phase delay with respect to the mag- 
netic field pulses. The time resolution is given by the 
x-ray pulse width and the jitter of the electronics, and 
is estimated as 100 ps. The domain wall motion in the 
two ferromagnetic (FM) layers was observed separately 
by tuning the photon energy of circularly polarized x- 
rays to the Fe and Co I/3 absorption edges for addressing 
the FeNi and Co layer, respectively. Measurements were 
performed at the UE56/2-PGM2 beamhne at BESSY in 
Berlin. 

In Fig. n we show images of the nucleation of domains 
and their subsequent expansion in the FeNi layer, ob- 
served at an amplitude of the bipolar field pulses of 5.8 
mT. The pulse shape is shown in Fig. ^ (e). The pos- 
itive and negative parts have the same amplitude (5.8 
mT) and length (47 ns). Figs. [H (a), (b), (c), and (d) 
show the magnetic domain structures in the FeNi layer 
at different pump-probe delay times, 9, 19, 24 and 29 
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FIG. 1: Bipolar magnetic field pulses (e) were used to reverse 
the magnetization in part of the FeNi layer. The amplitude 
of the field pulses was 5.80 mT on the plateaus of both the 
positive and negative part of the pulses. The magnetic do- 
main structure in the FeNi layer at different delay times is 
shown in (a-d), where the photon energy of circularly polar- 
ized x-rays was tuned to the Fe L3 absorption edge. White 
domains in (b) to (d) are discretized and superimposed to (f ) . 
Different shades of grey in (f) indicate different delay times. 
Images taken at the C0-L3 edge confirmed that no magneti- 
zation changes took place in the Co layer. 



ns, respectively. The easy axis is parallel to the vertical 
direction in the images, and the black and white con- 
trasts indicate the magnetization pointing down and up, 
respectively. The direction of the incoming x-rays with 
30° grazing angle to the surface is also indicated in (d). 
The field of view is 120 /xm. 

Before applying the pulses, the two FM layers were 
magnetically saturated in the same direction by a 15 mT 
field pulse of 3 ms duration. The direction is indicated by 
the white arrow in Fig. ^(d), and corresponds to black 
contrast. During the pump-probe experiments, the first, 
negative part of the pulse is in the same direction as 
the saturation, and was used to obtain the same starting 
condition for each pulse. The second, positive part is in 
the opposite direction (see (e)). The magnetization in 
the FeNi layer starts to reverse into the white direction 
during this second part of the pulse. The start of the 
positive part of the pulse is set to zero time. The plateau 
of the pulse starts at 9 ns, at which time one can see 
only some small fuzzy white splotches, probably due to 
magnetic domains that are blocked on surface defects (a). 
However, they do not expand with time, and are thus not 
regarded in the analysis. Some white domains become 
visible in (b), and, as time goes by, they get larger (c 
and d). To get sufficient statistics, millions of pulses were 
applied to acquire each image. 

The areas of white domains in Fig.Q] (b), (c), and (d) 
are superimposed into (f) with different shades of grey, 
middle grey, light grey, and dark grey, respectively. The 
expansion of the domains by propagation of the domain 
walls is clearly seen. In this experiment, the magnetiza- 



3 




Time (ns) 



FIG. 2: Increase of the domain perimeter vs. time during 
the plateau of the positive part of the field pulses for several 
field amplitudes. Broken lines are linear fits to the data. A 
simulation of the perimeter extension with time is shown by 
solid curves. 
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FIG. 3: Speed of perimeter extension vs. noHext (bottom 
axis) and HoiHext — Hcopi) (top axis). The linear fit to six of 
the data points ranging from 4.4 mT to 5.2 mT of fioiHext — 
Hcopi) indicates that the extension of the perimeter (i.e., the 
wall motion) is in the viscous regime. The exponential part 
at lower fields {fioHtot < 4.4 mT) is drawn using Eq. |5| In 
the higher field region, above Hwaiker, the perimeter speed 
is assumed to be constant (9.3 /im/ns). 



tion of the Co layer was not affected by the field pulses, 
i.e., it always showed the black contrast. 

In Fig. 13 the average perimeter of the three domains 
indicated by arrows in Fig. ^ (f ) is plotted as a fmiction 
of time for each amplitudes Hpuise of the field. ^oHext is 
ranging from 5.62 to 6.81 mT. The perimeter P is mea- 
sured by an ellipse-fit to the domains with an error of 
about 10%. The three domains were nucleated for all am- 
plitudes of the field pulses, and it was possible to measure 
P before they merged with other surrounding domains. 
The start of the plateau is indicated by a vertical line 
in Fig. 121 Linear fits to the data (broken lines) indicate 
that the measured values for P extended linearly with 
time for each field amplitude. However, the intersections 
to the time axis of these lines do not come to the same 
point, but show an apparent delay in the appearance of 
the domains, which amounts up to 17 ns for H),xt = 5.62 
mT. 

In Fig. 131 the slope of linear fits to the experimental 
data of Fig. (21 i.e., the speed of perimeter expansion for 
large domains, is plotted vs. iioHext (bottom axis) and 
vs. Ho{Hext- Hcopi) (top axis). In the range of fJ-oHext be- 
tween 5.6 and 6.4 mT, dP/dt increases linearly with field, 
meaning that the wall motion is in the viscous regime in 
this field range. A mobility of perimeter extension, fj,, of 
7800 m/(s mT) ±10% was obtained from the slope of a 
linear fit. The two data points with the highest ampli- 
tudes of the field pulses (6.70 mT and 6.81 mT) are not 
included in the linear fit, since the velocity seems to reach 
saturation at around 5.2 mT of fJ.o{Hext — Hcopi)- This 
saturation of velocity is due to the modification of the 
spin structure inside the domain wall above the so-called 
Walker limit field ^E^. 

For the analysis according to Eq. we use the data of 
Fig. 121 to extract the dependence of dP/dt on Htot, where 



we assume Htot ~ Hext — Hcopi for large domains. The 
negative contribution of H^aii to Htot for small domains 
makes it necessary to extend the velocity data to the 
lower field range for the numerical evaluation of P{t). 
We assumed thermally activated wall motion for lower 
Htot with perimeter extension velocity 



vp = vpQ exp{ 



fJ-oMsHtotVB — Ej^ 



(2) 



where Vb and Ep are the Barkhausen volume and the 
barrier energy, respectively [llli3,and /ioAfs = lTfor 
FeNi. To provide a smooth transition to the experimen- 
tally determined regime of viscous motion at H* = 4.4 
mT, values of Vb = (1.2 ± 0.1) x lO^^s ni^ and Ep = 
(2.6 ± 0.4) X 10-20 J 0] have to be used (solid line in 
Fig. 3), while wp,o of 55 nm/ns At Htot > 5.2 mT, 
a constant perimeter expansion speed of 9.3 /xm/ns was 
assumed in order to match the two data points at 6.70 
and 6.81 mT. 

The perimeter extension was numerically evaluated ac- 
cording to Eq. n with 1 ps time increments [l^ and 
is plotted in Fig. by solid curves. The best agree- 
ment was obtained for P{t = 0) = 1.9 /ini and 7A = 9 
mJ/m^. The simulation fits quite well to the experimen- 
tal data. P{t = 0) and 7 A are found to be correlated by 
7A = 3 X 10^ (mJ/m^) P{t ^ 0) + 3.2 (mJ/m^), and no 
significant worsening of the fit is observed for values of 
7A between 7 and 11 mJ/m^, if P{t = 0) is changed ac- 
cordingly and kept within 3% of the respective value. It 
is seen that a deviation from the linear behavior occurs 
at the initial stages of domain expansion, thus leading 
to the observed apparent delay. For higher ^^Hext (6.70 
mT and 6.81 mT), P starts to expand already on the 
rising slope of the pulse before 9 ns, as was also observed 
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in the experiment. 

The above simulation assumed that there is no real 
delay in domain nucleation. Since the grey scale con- 
trasts obtained in a static experiment were identical to 
the ones of the images obtained in the above experiment, 
the reproducibility of domain nucleation is close to 100%. 
Once the domains were visible, we did not observe an in- 
crease of the number of domains during the field pulses. 
This indicates that the domains were always nucleated as 
soon as the pulse field reached a certain threshold. The 
observation of a reproducible nucleation of domains at 
certain locations of the sample means that topographic 
features act as nucleation centers, and that the nucle- 
ation process is not thermally activated. Reproducible 
nucleation takes place at those positions where the local 
nucleation barrier is lower than the effective field. The 
increase of the number of nucleated domains with in- 
creasing pulse height indicates therefore a distribution of 
barrier heights. 

We have assumed a constant elliptical shape for the do- 
mains with aspect ratio 1:4, thus A « 24. The perimeter 
speed can be converted into a domain wall speed along 
the easy axis of magnetization. This is plotted on the 
right axis of Fig. [3 The mobility of viscous domain wall 
motion along the easy axis is then 1 600 m/ (s mT) ±10%. 
The initial P of 1.9 fj,m means that the width of the nu- 
cleation centers is around 0.19 /xm, which is too small to 
be seen in the images, the lateral resolution of which was 
set to w 1 ^m for these measurements. The maximum 
speed of 2000 m/s along the easy axis is an order of mag- 
nitude larger that the value found for narrow permalloy 
lines using similar field values, indicating that it might be 
possible to increase this speed by an appropriate shaping 
of the wires |0| . 

The Walker limit field is given by HoHi^aiker — 
l/2afioMs [13. For Fe2oNi8o, a value a = 0.01 is often 
used in the literature for the damping constant |3. l20ll2l| . 
With this value, HoHwaiker would be 5 mT. This is in 
good agreement with our experimentally obtained value 
of 5.2 mT. 

The constant effective domain wall energy 7 — 0.38 
mJ/ni^ that was used for the simulation includes ex- 
change and anisotropy energy inside the wall, magne- 
tostatic contributions from the sharp ends of the elliptic 
domains, as well as repulsion of the long sections of the 
magnetic domain walls 22] . A refined model should take 
into account the dependence of 7 on P. Since in the sim- 
ulation 7 contributes mainly when domains are small, 
deviations in the shape of the domains may also play a 
role when comparing the value for 7 to an estimate for 
the wall energy of an isolated Neel type wall, which leads 
to a much higher value of 2.3 mJ/m^ [23| . 

In conclusion, the reproducible nucleation of domains 
in the soft layer of a tunnel junction-like trilayer and the 
subsequent domain growth by domain wall propagation 
in the viscous regime and up to the Walker limit field is 



significantly influenced by the energy of magnetic domain 
walls. This leads to an apparent delay of domain expan- 
sion. A faster magnetization reversal, which is required 
for the coming magnetic recording devices, is drastically 
hindered by this magnetic domain wall energy when mag- 
netic domains are small. By analogy, the techniques and 
results shown in this work are equally applicable to other 
nucleation and growth phenomena like crystallographic 
phase changes, superconducting transitions, and solidifi- 
cation theory. 
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